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Abstract. We study the spatial patterns formed by inter- 
acting populations or reacting chemicals under the influence 
of chaotic flows. In particular, we have considered a three- 
component model of plankton dynamics advected by a mean- 
dering jet. We report general results, stressing the existence 
of a smooth-filamental transition in the concentration pat- 
terns depending on the relative strength of the stirring by the 
chaotic flow and the relaxation properties of planktonic dy- 
namical system. Patterns obtained in open and closed flows 
are compared. 



1 Introduction 

The transport of biologically or chemically active substances 
by a fluid flow is a problem of great geophysical relevance. 
Important examples arise in the study of atmospheric advec- 
tion of reactive pollutants or chemicals, such as ozone, N2O 
(Mahlman et al., 1984), or in the dynamics of plankton pop- 
ulations in ocean currents (Denman and Gargett, 1995). The 
inhomogeneous nature of the resulting spatial distributions 
was recognized some time ago (Steele, 1978, and references 
therein). More recently, satellite remote sensing and detailed 
numerical simulations identify filaments, irregular patches, 
sharp gradients, and other complex structures involving a 
wide range of spatial scales in the concentration patterns. 
In the case of atmospheric chemistry, the presence of strong 
concentration gradients has been shown to have profound im- 
pact on global chemical time-scales (Edouard et al., 1996). 
On-site measurements and data analysis of the chemical or 
biological fields have confirmed their fractal or multifractal 
character (Pascual et al., 1995; Seuront et al., 1996; Bacmeis- 
ter et al., 1997; Tuck and Hovde, 1999). 

In the case of plankton communities, patchiness has been 
variously attributed to the interplay of diffusion and biolog- 
ical growth, oceanic turbulence, diffusive instabilities, and 
nutrient or biological inhomogeneities (Mackas et al., 1985). 
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Advection by unsteady fluid flow and predator-prey interac- 
tions (formally equivalent to chemical reaction) are emerging 
as two key ingredients able to reproduce the main qualitative 
features of plankton patchiness (Abraham, 1998). 

The 'chaotic advection' paradigm has been shown to be 
a useful approach to understand geophysical transport pro- 
cesses at large scales (Haynes, 1999). Briefly, chaotic ad- 
vection (or Lagrangian chaos)(Aref, 1984) refers to the La- 
grangian complex motion of fluid parcels arising from a flow 
which is not turbulent in the Eulerian description. Lagrangian 
chaotic flows are much simpler than turbulent ones, being 
thus more accessible to analytical characterization and un- 
derstanding. They retain however many of the qualitative 
features relevant to transport and mixing processes in com- 
plex geophysical flows. 

Though the properties of inert passive tracer fields under 
chaotic advection have been widely studied (Ott and Anton- 
sen, 1989) much less is known about biologically or chemi- 
cally evolving reactant distributions. Nonetheless, some re- 
sults have been recently obtained, as for example in reactions 
of the type A + B —* C in closed chaotic flows (Metcalfe 
and Ottino, 1994) and in open chaotic flows (Toroczkai et al., 
1998). Recently, some of us (Neufeld et al., 1999a) consid- 
ered the general case of stable chemical dynamics in closed 
chaotic flows in the limit of small diffusion and in the pres- 
ence of an external spatially non-homogeneous source of one 
of the chemical components. The main result was that the re- 
lationship between the rate at which the chemical dynamics 
approaches local equilibrium with the chemical source and 
the characteristic time scale of the stirring by the chaotic flow 
determines the fractal or non-fractal character of the long- 
time distribution. The faster the stirring is, the more irregular 
is the pattern. 

The purpose of this Paper is to apply and verify the gen- 
eral results above in a concrete model of plankton dynamics 
in flows of geophysical relevance. In addition we compare 
structures appearing in closed and open flows, stressing the 
intermittent character of irregularities in the open flow case. 
We expect this result to apply also to other situations in at- 



mospheric or oceanic chemistry. 

The paper is organized as follows: next Section summa- 
rizes the general theoretical results obtained by Neufeld et al. 
(1999a). The particular plankton dynamics and the two dif- 
ferent flows subject of our study are presented in Sect. 3. 
They are variations of a kinematic model for a two-dimensional 
meandering jet, leading to a closed and an open flow model. 
Section 4 describes numerical results for the closed flow case, 
whereas Sect. 5 considers the open flow. Finally, Sect. 6 con- 
tains our conclusions. 



2 General results 

The temporal evolution of reacting fields is determined by 
advection-reaction-diffusion equations. Advection because 
they are under the influence of a flow, reaction because we 
consider species interacting with themselves and/or with the 
carrying medium. Diffusion because turbulent or molecular 
random motion smoothes out the smallest scales. For the 
case of an incompressible velocity field v(r, t), the standard 
form of these equations is 



dd(T,t) 
dt 



v(r,t)-VCi(r,t) = 

MC 1 ,...,C Nl r) + vV 2 C t (r,t) (1) 



where C,(r, t), i = 1, ...N, are interacting chemical or bio- 
logical fields advected by the flow v(r, t), /,(Ci, ..., Cn, r) 
are the functions accounting for the interaction of the fields 
(e.g. chemical reactions or predator-prey interactions). Dif- 
fusion effects are only important at small scales and we will 
neglect them in the following. In this limit of zero diffusion 
v — ► the above description can be recast in Lagrangian 
form: 






da 

dt 



= fi(C u C 2 ,..,C N ,v = r(t)), i = l,.,N, 



(2) 



(3) 



where the second set of equations describes the chemical or 
population dynamics inside a fluid parcel that is being ad- 
vected by the flow according to the first equation. In the 
absence of diffusion, a coupling between the flow and the 
chemical/biological evolution can only appear as a conse- 
quence of the spatial dependence of the fi(C\, ..., Cn,t) 
functions. This spatial dependence describes non-homogeneous 
sources or sinks for the chemical reactants or spatially non- 
homogeneous reaction or reproduction rates. Such inhomo- 
geneities may arise naturally from a variety of processes such 
as localized upwelling, inhomogeneous solar irradiation, or 
river run-off, to name a few. 

From now on, the incompressible flow v(r, t) will be as- 
sumed to be two-dimensional and time dependent. This situ- 
ation generally leads to chaotic advection. For simplicity, our 
general arguments will be stated for the case in which v(r, t) 
satisfies the technical requirement of hyperbolicity, but in the 



examples less restrictive flows will be used. The most salient 
feature of advection by a chaotic flow is sensibility to ini- 
tial conditions, that is, fluid particles initially close typically 
diverge in time at a rate given by the maximum Lyapunov 
exponent of the flow Xf > 0: 



\6r(t)\ ~ |<5r(0)|e 



x F t 



(4) 



Equation (4) is valid for nearly all the initial orientations of 
the initial particle separation <5r(0). However, the incom- 
pressibility condition implies that there is a particular orien- 
tation of the initial separations <5r c (0) for which the two tra- 
jectories approach each other: |<5r(£)| ~ |<5r c (0)|e Aj7 ', with 
A F = — Xf- 

The general class of chemical reactions studied by Neufeld 
et al. (1999a) was the one leading to stable local equilibrium 
in the absence of flow or, in terms of the chemical dynami- 
cal subsystem (3), the dynamics approaching a unique fixed 
point for each constant position r. This means that there are 
not chemical instabilities nor chemical chaos, and that the 
concentrations tend to approach a value determined at each 
point by the sources in /j. In the presence of advection by 
the flow r = f (£), the relaxation process is altered, but can 
be characterized by the value of the maximum Lyapunov ex- 
ponent Ac of the chemical subsystem (3), which we assume 
to remain negative. 

For characterizing the spatial structure of the Cj fields we 
calculate the difference 



8C i = C i (r + 6r,t)-C i (r,t) 



(5) 



Inserting this expression, for \6r\ small enough, into the equa- 
tion for the chemical dynamics (details can be found in Neufeld 
et al., 1999a) we can obtain the evolution of the gradients of 
the chemical field at long times: 



VCi(r,t) 



N 

E 



v(c? • v?)v?e< A * +Ao >' 



ds Ef=i {vj 3 (v( S ))-^m^ XF+Xc){t - 



.(6) 



where the vectors V* = n (f (£)) m are combinations of the 
vectors n(r) pointing at each point along the most contract- 
ing direction of the flow and of the vector {vi,i = 1, ..., N} 
associated to the contracting direction in the chemical sub- 
space, f (s) is the trajectory ending at r at time t. It is im- 
portant to realize that Eq. (6) gives the long time behavior 
of VCi correctly in all but in one direction. The directional 
derivative of C\ in the most expanding direction should be 
obtained with (6) but replacing Xp by X' F = —Xp, and the 
vectors n by the ones associated to the expanding direction. 
The convergence of the gradients for t — > oo depends on 
the sign of the exponent Xp+Xc- There are two possibilities: 

- If Xp + Xc < then the convergence of the chemical 
dynamics towards local equilibrium is stronger than the 
effect of the chaotic flow on the fluid particles. Gradi- 
ents are finite so that a smooth asymptotic distribution 
is attained by the chemical fields. 



If Xp + Xc > then in the t — > oo limit the chemical 
pattern becomes nowhere differentiable. An irregular 
structure with fractal properties is developed. Remem- 
ber however that at each point there is a direction for 
which Xp + Xc should be substituted by — Xp + Xc, 
always negative, in (6). In this direction derivatives are 
finite and the field is smooth. 

Thus the resulting structure is Mamental, i.e., irregular 
in all directions except in one along which it is smooth. 
This one corresponds to the direction of the filaments ly- 
ing along the unstable foliation of the chaotic advection. 
The fractal characteristics of the filamental structure in 
the closed flow case were investigated by Neufeld et al. 
(1999a). We will see however that there are differences 
between the closed and the open case, which will be 
discussed in Sects. 4 and 5. 



3 The plankton and the jet models 

In the numerical investigations below we will consider a sim- 
ple model of plankton dynamics immersed in a meandering 
jet flow. 

This plankton model, used by Abraham (1998) and related 
to the one used by Levin and Segel (1976), considers ex- 
plicitly three trophic levels: the nutrient content of a water 
parcel, described in terms of its carrying capacity C (defined 
as the maximum phytoplankton content it can support in the 
absence of grazing), the phytoplankton biomass P, and the 
zooplankton Z. The Lagrangian 'chemical' subsystem (3) 
reads: 



1 — tanh 



dC 

~dt 
dP 

~dt 
dZ 
~dt 



a(C (r)-C) 
PZ - 5Z 2 . 



(7) 
(8) 
(9) 



All terms have been adimensionalized to keep a minimal 
number of parameters. Equation (7) describes the relaxation 
of the carrying capacity, at a rate a, towards an inhomoge- 
neous shape Cb(r). This will be the only explicitly inhomo- 
geneous term in the model, and describes a spatially depen- 
dent nutrient input, arising from some topography-determined 
upwelling distribution or latitude dependent illumination, for 
example. The first terms in Eq. (8) describe phytoplankton 
logistic growth, whereas the last one models predation by 
zooplankton. This effect gives also rise to the first term in 
(9). The term containing S, the zooplankton mortality, de- 
scribes zooplankton death produced by higher trophic lev- 
els. The only stable fixed point of model (7)-(9) is given by 
C* = C (r), P* = C 5/(5 + C ), and Z* = P*/S. 

The model flow will be given by the following streamfunc- 
tion (Bower, 1991): 

<P(x, y) = 



y — B(t) cos [k(x — ct)] 
(1 + k 2 B(t) 2 sin 2 [k(x - ct)])- 



(10) 



It describes a jet flowing eastwards, with meanders in the 
North-South direction which are themselves advected by the 
jet at a phase velocity c. B(t) and k are the (properly adimen- 
sionalized) amplitude and wavenumber of the undulation in 
the streamfunction. 

The motion of the tracer particles (the dynamical system 
(2) ) is given by 



(11) 



If the amplitude B of the meanders is time-independent, 
a simple change in the frame of reference renders the flow 
time-independent and Eq. (11) defines a non-chaotic inte- 
grable dynamical system. Chaotic advection appears in this 
model if B is made to vary in time, for example periodically: 



dx 


d(j> 


dt 


dy 


dy 


d<j> 


dt 


dx 



B(t) = B + ecos(tot 



(12) 



Following Cencini et al. (1998), we use the parameter val- 
ues B = 1.2, c = 0.12, k = 2n/L x ,L x = 7.5, w = 
0.4, e = 0.3 and 6 = \. These values guarantee the ex- 
istence of 'large scale chaos', i.e, the possibility that a test 
particle crosses the jet passing from North to South or vicev- 
ersa. This is weaker than the requirement of hyperbolicity, 
but is enough to illustrate the general aspects of our theory. 

The natural interpretation of the jet-flow just introduced is 
as an open flow: it advects most of the fluid particles from 
x = — oo towards x = oo. We will localize the source term 
near the origin of coordinates: 



Oo{x,y) 



1 + Asm(2irx/ L x ) siii(2iry / L y ) 
if — L x < x < L x 

elsewhere 



(13) 



In this way, there is inhomogeneous nutrient input just near 
the origin, and capacity and plankton concentration in the 
fluid particles will decay as they are advected downstream. 

A quite different class of natural flows are closed ones, e.g. 
recirculating flows in closed basins. Our jet model flow can 
be made closed simply by imposing periodic boundary con- 
ditions at the ends of the interval — L x < x < L x . Particles 
leaving the region through the right boundary are reinjected 
from the left. In this way nutrients are injected and extracted 
continuously from fluid elements as they traverse the differ- 
ent regions of the source 



Co{x,y) = 1 + Asui(2ttx/ L x ) sin(2TT y/Ly) 



(14) 



This was the situation considered by Neufeld et al. (1999a). 
We will see that different structures develop in the open and 
in the closed situation. 
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Fig. 1. Phytoplankton smooth pattern in the closed flow, with an horizontal 
transect taken along y = 0.8 



Fig. 2. Phytoplankton filamental pattern in the closed flow, with a transect 
along y = 0.8 



4 Closed flows 



Numerically we proceed by integrating backwards in time 
Eq. (2) with initial coordinates on a rectangular grid (300 x 
150) and then the chemical field for each point is obtained by 
integrating (3) forward in time along the fluid trajectories so 
generated. 

Figure 1 shows a snapshot of the long-time phytoplankton 
distribution in the closed flow case for parameter values a = 
0.25, S = 2.0, L y = 4. and A = 0.2. In this case Xp + X c < 
0, therefore the distribution is smooth. A transect along the 
line y — 0.8 is also shown. Taking a = 0.025, so that now 
Xf + Xc > we obtain the distribution in Fig. 2. A complex 
filamental structure is clearly seen, in agreement with our 
theoretical arguments. The fractal nature of the pattern is 
also seen in the horizontal cut presented also in Fig. 2. A 
Holder exponent of Xc\/Xp was predicted for these kind of 
transects in (Neufeld et al., 1999a). This implies plankton- 
variance power spectrum decaying as fc _/3 , with /3 = 1 + 
2\Xc\/Xf- Thus j3 is in the range ]1, 3] which agrees with 
field observations of plankton distributions (Seuront et al., 
1996). 



5 Open flows 

Contrarily to closed flows, open flows are characterized by 
unbounded trajectories of fluid particles. Typical fluid parti- 
cles enter and, after some time, leave the region of the nu- 
trient source (Co(r) > 0). Thus, most fluid elements have 
only spent the most recent part of their trajectories inside the 
active region, with the rest of their evolution spent in regions 
with no spatial dependence of the nutrient input. For this part 
of the trajectory the gradient V/j in Eq. (6) vanishes, imply- 
ing that there would not be divergence of the gradients in the 
long-time limit. 

It is well known from the study of chaotic advection in 
open flows (Pentek et al., 1995), that while most of the par- 
ticles spend only a finite amount of time in selected bounded 



regions of the flow, typically these regions contain also bounded 
orbits in which some of the particles can stay forever. Al- 
though the chaotic set formed by the bounded orbits is a frac- 
tal set of measure zero, particles visiting the vicinity of the 
stable manifold of this set can still spend arbitrarily long time 
in the selected bounded region. This leads to the formation of 
characteristic fractal patterns in the advection dynamics even 
in the case of passive particles, as it was shown in numeri- 
cal studies (Jung et al., 1993) and laboratory experiments of 
open flows (Sommerer et al., 1996). 

For the particles that have spent infinitely long time in the 
source region, their chemical/biological evolution is equiv- 
alent to the one in a closed flow with the possibility of di- 
verging gradients. The only difference is that now the val- 
ues of Xp and Xc are those corresponding to the chaotic set 
of bounded orbits that never leave the biologically active re- 
gion. The smooth-filamental transition observed in the pre- 
vious Section for the closed flow will occur here only on this 
fractal set, which will always be surrounded by a smooth dis- 
tribution. This results in a strongly intermittent character of 
the filamental field. Figure 3 shows a phytoplankton pattern 
for parameter values a = 0.025, <5 = 2.0 and A = 0.2, 
and a transect crossing it, in the open-flow case. The inho- 
mogeneity in the filamental structure is obvious, with singu- 
larities in a set recognized as the unstable manifold of the 
chaotic set formed by bounded orbits in the source region. 
Smooth structures are also obtained when Xp + Xc < 0. 
The proper description of the resulting structures should use 
the concept of multifractality, that is, inhomogeneous distri- 
bution of fractal properties. In fact, even in the closed flow 
case, at finite times the flow Lyapunov exponent Xp would 
have space-dependent finite-time corrections, that will need 
to be taken into account in a proper generalization of (6). A 
quantitative description of these multifractal filamental struc- 
tures will be presented elsewhere (Neufeld et al., 1999b). 
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Fig. 3. Phytoplankton filamental pattern in the open flow, with a transect 
along y = 0.8. 



6 Conclusions 

The spatial patterns formed by interacting populations under 
the influence of chaotic flows have been studied. In particu- 
lar, we have considered a coupled model of 'nutrient', phy- 
toplankton and zooplankton concentrations advected by two 
(open and closed) jet-like flows. General results have been 
reported for arbitrary chaotic flows, stressing the existence 
of a smooth-filamental transition depending on the relative 
strength of the maximum Ly apunov exponent of the flow and 
the one corresponding to the planktonic dynamical system. 
Patterns obtained for open and closed chaotic flows are dif- 
ferent because of the transient or permanent character of the 
biological activity. Comparison of the structures for the dif- 
ferent fields (capacity, phytoplankton and zooplankton), and 
quantitative description of their multifractal properties will 
be presented elsewhere (Neufeld et al., 1999b). 

The models considered here are extreme simplifications of 
real biological and geophysical situations. We expect how- 
ever that the main qualitative features found here, namely 
the possibility of finding smooth or filamental patterns de- 
pending on stirring and relaxation rates, and the increased 
inhomogeneities in open flows, to be present in more realis- 
tic chemical or biological transport situations. 
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